Parallel Batch-Dynamic Graph Connectivity

In this paper, we study batch parallel algorithms for the dynamic connectivity problem, a fundamental problem that has received considerable attention in the sequential setting. The most well known sequential algorithm for dynamic connectivity is the elegant level-set algorithm of Holm, de Lichtenberg and Thorup (HDT), which achieves $O(\log^2 n)$ amortized time per edge insertion or deletion, and $O(\log n / \log\log n)$ time per query. We design a parallel batch-dynamic connectivity algorithm that is work-efficient with respect to the HDT algorithm for small batch sizes, and is asymptotically faster when the average batch size is sufficiently large. Given a sequence of batched updates, where $\Delta$ is the average batch size of all deletions, our algorithm achieves $O(\log n \log(1 + n / \Delta))$ expected amortized work per edge insertion and deletion and $O(\log^3 n)$ depth w.h.p. Our algorithm answers a batch of $k$ connectivity queries in $O(k \log(1 + n/k))$ expected work and $O(\log n)$ depth w.h.p. To the best of our knowledge, our algorithm is the first parallel batch-dynamic algorithm for connectivity.Comment: This is the full version of the paper appearing in the ACM Symposium on Parallelism in Algorithms and Architectures (SPAA), 201

BLITZEN: A highly integrated massively parallel machine

The architecture and VLSI design of a new massively parallel processing array chip are described. The BLITZEN processing element array chip, which contains 1.1 million transistors, serves as the basis for a highly integrated, miniaturized, high-performance, massively parallel machine that is currently under development. Each processing element has 1K bits of static RAM and performs bit-serial processing with functional elements for arithmetic, logic, and shifting

A Paradox of State-Dependent Diffusion and How to Resolve It

Consider a particle diffusing in a confined volume which is divided into two equal regions. In one region the diffusion coefficient is twice the value of the diffusion coefficient in the other region. Will the particle spend equal proportions of time in the two regions in the long term? Statistical mechanics would suggest yes, since the number of accessible states in each region is presumably the same. However, another line of reasoning suggests that the particle should spend less time in the region with faster diffusion, since it will exit that region more quickly. We demonstrate with a simple microscopic model system that both predictions are consistent with the information given. Thus, specifying the diffusion rate as a function of position is not enough to characterize the behaviour of a system, even assuming the absence of external forces. We propose an alternative framework for modelling diffusive dynamics in which both the diffusion rate and equilibrium probability density for the position of the particle are specified by the modeller. We introduce a numerical method for simulating dynamics in our framework that samples from the equilibrium probability density exactly and is suitable for discontinuous diffusion coefficients.Comment: 21 pages, 6 figures. Second round of revisions. This is the version that will appear in Proc Roy So

High-resolution investigations of ripple structures formed by femtosecond laser irradiation of silicon

We report on the structural investigation of self-organized periodic microstructures (ripples) generated in Si(100) targets after multishot irradiation by approximately 100-fs to 800-nm laser pulses at intensities near the single shot ablation threshold. Inspection by surface sensitive microscopy, e.g., atomic force microscopy (AFM) or scanning electron microscopy (SEM), and conventional and high-resolution transmission electron microscopy reveal complex structural modifications upon interaction with the laser: even well outside the ablated area, the target surface exhibits fine ripple-like undulations, consisting of alternating crystalline and amorphous silicon. Inside the heavily modified area, amorphous silicon is found only in the valleys but not on the crests which, instead, consist of highly distorted crystalline phases, rich in defects

Self-replication and evolution of DNA crystals

Is it possible to create a simple physical system that is capable of replicating itself? Can such a system evolve interesting behaviors, thus allowing it to adapt to a wide range of environments? This paper presents a design for such a replicator constructed exclusively from synthetic DNA. The basis for the replicator is crystal growth: information is stored in the spatial arrangement of monomers and copied from layer to layer by templating. Replication is achieved by fragmentation of crystals, which produces new crystals that carry the same information. Crystal replication avoids intrinsic problems associated with template-directed mechanisms for replication of one-dimensional polymers. A key innovation of our work is that by using programmable DNA tiles as the crystal monomers, we can design crystal growth processes that apply interesting selective pressures to the evolving sequences. While evolution requires that copying occur with high accuracy, we show how to adapt error-correction techniques from algorithmic self-assembly to lower the replication error rate as much as is required

Thermodynamic entropy of a many body energy eigenstate

It is argued that a typical many body energy eigenstate has a well defined thermodynamic entropy and that individual eigenstates possess thermodynamic characteristics analogous to those of generic isolated systems. We examine large systems with eigenstate energies equivalent to finite temperatures. When quasi-static evolution of a system is adiabatic (in the quantum mechanical sense), two coupled subsystems can transfer heat from one subsystem to another yet remain in an energy eigenstate. To explicitly construct the entropy from the wave function, degrees of freedom are divided into two unequal parts. It is argued that the entanglement entropy between these two subsystems is the thermodynamic entropy per degree of freedom for the smaller subsystem. This is done by tracing over the larger subsystem to obtain a density matrix, and calculating the diagonal and off-diagonal contributions to the entanglement entropy.Comment: 18 page

Electronic Theory for the Nonlinear Magneto-Optical Response of Transition-Metals at Surfaces and Interfaces: Dependence of the Kerr-Rotation on Polarization and on the Magnetic Easy Axis

We extend our previous study of the polarization dependence of the nonlinear optical response to the case of magnetic surfaces and buried magnetic interfaces. We calculate for the longitudinal and polar configuration the nonlinear magneto-optical Kerr rotation angle. In particular, we show which tensor elements of the susceptibilities are involved in the enhancement of the Kerr rotation in nonlinear optics for different configurations and we demonstrate by a detailed analysis how the direction of the magnetization and thus the easy axis at surfaces and buried interfaces can be determined from the polarization dependence of the nonlinear magneto-optical response, since the nonlinear Kerr rotation is sensitive to the electromagnetic field components instead of merely the intensities. We also prove from the microscopic treatment of spin-orbit coupling that there is an intrinsic phase difference of 90$^{\circ }$ between tensor elements which are even or odd under magnetization reversal in contrast to linear magneto-optics. Finally, we compare our results with several experiments on Co/Cu films and on Co/Au and Fe/Cr multilayers. We conclude that the nonlinear magneto-optical Kerr-effect determines uniquely the magnetic structure and in particular the magnetic easy axis in films and at multilayer interfaces.Comment: 23 pages Revtex, preprintstyle, 2 uuencoded figure

Restricted random walk model as a new testing ground for the applicability of q-statistics

We present exact results obtained from Master Equations for the probability function P(y,T) of sums $y=\sum_{t=1}^T x_t$ of the positions x_t of a discrete random walker restricted to the set of integers between -L and L. We study the asymptotic properties for large values of L and T. For a set of position dependent transition probabilities the functional form of P(y,T) is with very high precision represented by q-Gaussians when T assumes a certain value $T^*\propto L^2$. The domain of y values for which the q-Gaussian apply diverges with L. The fit to a q-Gaussian remains of very high quality even when the exponent $a$ of the transition probability g(x)=|x/L|^a+p with 0<p<<1 is different from 1, all though weak, but essential, deviation from the q-Gaussian does occur for $a\neq1$. To assess the role of correlations we compare the T dependence of P(y,T) for the restricted random walker case with the equivalent dependence for a sum y of uncorrelated variables x each distributed according to 1/g(x).Comment: 5 pages, 7 figs, EPL (2011), in pres

Quantum heat engines and nonequilibrium temperature

A pair of two-level systems initially prepared in different thermal states and coupled to an external reversible work source, do not in general reach a common temperature at the end of a unitary work extraction process. We define an effective temperature for the final nonequilibrium but passive state of the bipartite quantum system and analyse its properties.Comment: Five pages, Accepted for publication in Physical Review

Landau-Drude Diamagnetism: Fluctuation, Dissipation and Decoherence

Starting from a quantum Langevin equation (QLE) of a charged particle coupled to a heat bath in the presence of an external magnetic field, we present a fully dynamical calculation of the susceptibility tensor. We further evaluate the position autocorrelation function by using the Gibbs ensemble approach. This quantity is shown to be related to the imaginary part of the dynamical susceptibility, thereby validating the fluctuation-dissipation theorem in the context of dissipative diamagnetism. Finally we present an overview of coherence-to-decoherence transition in the realm of dissipative diamagnetism at zero temperature. The analysis underscores the importance of the details of the relevant physical quantity, as far as coherence to decoherence transition is concerned.Comment: 8 pages and 5 figure